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ABSTRACT 


The  reliability  of  a  system  is  the  probability  that  the 
system  will  survive  or  complete  an  intended  mission  of 
certain  duration. 

Describing  all  possible  ways  that  a  system  can  survive  a 
mission  in  reliability  shorthand  gives  a  simple  approach  to 
reliability  computations.  Reliability  computation  for  a 
system  defined  by  shorthand  notation  is  greatly  dependent 
upon  the  convolution  problem. 

Assuming  constant  component  failure  rates,  this  paper 
presents  an  analytical  approach  and  a  computer  program  for 
computing  the  reliability  of  any  convolution  of  independent 
and  exponentially  distributed  random  variables.^ 

r  . 

••  j  .  jfQx-  i  .  .  *•<• 
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I.  INTRODUCTION 


In  making  a  mathematical  model  for  a  real-life  phenom¬ 
enon  it  is  always  necessary  to  make  certain  simplifying 
assumptions  so  as  to  render  the  mathematics  tractable.  One 
of  the  simplifying  assumptions  that  is  often  made  is  to 
assume  that  certain  random  variables  are  exponentially 
distributed.  The  reason  for  this  is  that  the  exponential 
distribution  is  relatively  easy  to  work  with. 

The  property  of  the  exponential  life  distribution  which 
makes  it  easy  to  analyze  is  its  (memoryless)  lack  of  deteri¬ 
oration  with  time.  By  this  we  mean  that  if  the  life  time  of 
an  item  is  exponentially  distributed,  then  an  item  which  has 
been  in  use  for  a  certain  amount  of  time  is  as  good  as  a  new 
item  in  regards  to  the  amount  of  time  remaining  until  it 
fails.  When  the  life  time  of  an  item  is  exponentially 
distributed,  the  failure  rate  function  for  the  item  is 
constant . 

Under  the  assumption  of  constant  component  failure 
rates,  it  is  possible  to  build  a  reliability  shorthand 
[Ref.  1]  for  any  system.  The  term  system  is  used  to  describe 
a  set  of  components  organized  to  perform  some  mission. 


Any  study  on  system  reliability  requires  a  description 
for  the  system's  life  and  a  derivation  of  the  system's 
survival  function.  The  reliability  shorthand  gives  a  simple 
and  easy  way  for  describing  a  system's  life,  but  it  is 
difficult  to  implement  computationally  since  it  involves 
considerable  complexity  in  handling  convolutions.  Here,  the 
term  convolution  refers  to  the  summation  of  independent 
random  variables  (lives). 

This  paper  presents  an  analytical  approach  for  obtaining 
a  general  equation  for  the  survival  function  of  any  convolu¬ 
tion  of  independent  and  exponentially  distributed  random 
variables . 

Section  2  deals  with  convolutions  in  detail  and  gives  a 
mathematical  derivation  for  a  general  equation  for  the 
survival  function  for  any  convolution  of  exponential  random 
variables,  using  Laplace  transforms  and  Theorem  of  Residues. 

Section  3  gives  a  mathematical  derivation  for  an  alter¬ 
native  formula  for  the  survival  function  of  any  convolution 
of  exponential  random  variables  by  computing  the  coeffi¬ 
cients  of  all  the  polynomials  that  accompany  the  exponential 
terms  in  the  survival  function  equation. 

Appendix  A  contains  a  computer  program  written  in 
Fortran  for  computing  the  reliability  of  any  convolution  of 
inependent  and  exponentially  distributed  random  variables. 
This  program  uses  the  general  equation  in  Section  2. 
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II.  CONVOLUTIONS 


Under  the  assumption  of  constant  component  failure 
rates,  this  section  presents  an  analytical  approach  for 
finding  a  general  equation  for  the  survival  function  of  any 
convolution  of  independent  and  exponentially  distributed 
random  variables. 

I 

A.  SURVIVAL  FUNCTION  BY  INTEGRATION 

Let  us  start  with  a  simple  system;  a  standby  system 
having  one  active  component  and  one  cold  spare  component. 
The  life  time  of  the  system  is 

T  =  T,  +  Tz  , 

where 

T,  ~  EXP(X|  )  , 
t2^  exp(X2)  , 

and 

Tj,  T^  are  independent. 

An  active  component  A  is  to  complete  a  mission  of  dura¬ 
tion  t,  while  a  spare  component  S  replaces  the  active  compo¬ 
nent  when  it  fails.  The  life  duration  of  the  active 
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Figure 


WORKING  COMPONENTS  ' - 1 

2.1  A  Single  Active  Component  With  One  Spare 


component  is  T  j  and  the  life  duration  for  the  spare  compo 
nent,  if  it  is  used,  is  T2  (see  Fig  2.1). 

The  system  survival  function  is  given  by 

?T  (t)  =  P  (  T  >  t) 


/ 


F  T)  ( t )  +  /  FTz(t-s).fT4(s) 


ds 


=  e 


-X,t  f  -X,S 


*  J  *  K 


ds  , 


where 


fT)(s)  =  -  dFT^(s)/ds 

■  X,rx>s  . 


By  integration,  the  survival  function  for  the  system  is 


Now,  suppose  a  system  consists  of  one  active  component 
with  constant  failure  rate  Xj  and  two  spares  having  constant 
failure  rates  \  ^  and  X3  respectively.  The  life  time  of  the 
system  is 


where 


and 


T  =  T,  +  T2  +  T3  , 

T,  ~  EXP(X)  )  , 

Tz  ~  EXP(XZ  )  , 

T3  ~  EXP  (  X3  )  , 

T | ,  T  2  i  T  3  are  independent . 


The  shorthand  notation  for  the  systemi is 

EXP ( X ,  )  +  EXP(Xa)  +  EXP(X3)  • 

The  survival  function  of  the  sustera  is  given  by 


Fy  (t)  = 


where 


F-j^(t-s) 
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=  "dF-q  +  T^^/ds 

_Xj_Xz_r  -\2s  -\,s 

'  x, — x.2.  L 

By  integration,  the  survival  function  for  the  system  is 
FT  (t)  =  Fjj  +  ( t ) 

-  XzX3  *  ^]tl  X,  X3  e“  t  kjkz  e"  t>0 
(XfXi)(XrX0  (XfXt)(X j~\z)  (XfX3)(Xi-X3) 

To  compute  the  survival  function  for  a  system  of  n 
lives , 

T  =  T,  +  .. .♦  Tn  , 

we  can  proceed  in  a  similar  manner.  However  the  integrations 
will  become  increasingly  complex,  particularly  if  the  tacit 
assumption  made  so  far  that  all  component  failure  rates  are 
different  is  abandoned.  Therefore,  we  will  use  another  math¬ 
ematical  way  to  compute  the  survival  function  without  direct 
integration.  This  will  be  done  with  Laplace  transforms  and 
the  Theorem  of  Residues  in  the  next  subsections. 


B.  CLASSIFICATION  OF  CONVOLUTIONS 


Consider  a  system  which  has  a  life  time 
T  =  T,  +  ...  +  Tn  , 

where  T(  , . . . ,Tn  are  independent  and  exponentially  distrib¬ 
uted  random  variables.  This  system  can  be  classified  as  one 
of  three  possible  convolution  cases  : 


1. 

First  Case 

The  system 

has 

n  dissimilar 

failure  rates, 

X;  * 

>  i  H  ,  i,j 

c 

H 

II 

2. 

Second  Case 

The  system 

has 

n  identical 

failure  rates, 

11 

»  i»j  —  1,  •• 

.  ,  n. 

3. 

Third  Case 

The  system  has  a  combination  of  the  first  case  and 

«  • 

the  second  case,  some  similar  and  some  dissimilar  failure 
rates . 

We  consider  the  third  case  to  be  of  particular 
interest,  since  it  better  represents  the  general  situation 
in  computing  the  reliability  of  any  convolution  of  exponen¬ 
tially  distributed  random  variables. 


C.  CHARACTERISTIC  FUNCTIONS  AND  LAPLACE  TRANSFORMS 


Consider  a  random  variable  T  which  has  a  distribution 
function  F  j(t)  and  a  density  function  f j (t).  The  charac¬ 
teristic  function  of  T  is  defined,  for  any  real  number  u, 
[Ref.  2]  by 


+  o© 

r  iut 

J  ' 

-  oo 


f  T (t)  dt  • 


If  T  is  non  negative  (T  0)  and  exponentially  distributed 
with  a  failure  rate  \  ,  then 


• f  J  (t )  dt 


Now,  if 


T  =  T,  +  ...  +  Tn  , 

where 

Tj  —  EXP(\j  ),  i  =  1 . . 

and  T(,...,TJ-|  are  independent,  then  the  characteristic 
function  of  the  sum  (convolution)  Tj+  ...  +  Tn  is 


This  implies  that  the  characteristic  function  for  the  sum 
T,  +  ...  +  1,-j  is  just  the  product  of  the  individual  charac¬ 


teristic  functions. 


When  dealing  with  random  variables  which  only  assume  non 
negative  values,  it  is  more  convenient  to  use  Laplace  trans¬ 
forms  rather  than  characteristic  functions.  The  Laplace 
transform  of  the  random  variable  having  distribution  Fy(t) 
is  defined  [Ref.  2]  by 


,  r  -st  i 

<pT(s)  =  E  e 


oO 


-  st 

e  dF  T ( t ) 


■/ 


- St 

e  f  t (t)  dt , 


where  f-p(fc)  is  the  density  function  of  T.  This  integral 
exists  for  a  complex  variable  s=a+bi,  where  a>0 . 

If  Tj  , .  .  . ,Tn  are  independent  random  variables,  then  the 
Laplace  transform  of  the  sum  (convolution)  T,  +...+  Tn  is 
given  by 

-sCT  +  -  +  Tn) 


< (s )  =  E  e 


=  E 
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It  is  important  to  note  that  the  Laplace  transform  uniquely 
determines  the  distribution.  That  is,  in  fact,  there  is  a 
one  to  one  correspondence  between  distribution  functions  and 
Laplace  transforms. 

Now,  we  will  consider  the  following  question.  Given  a 
function  <^y(s),  does  there  exist  ,a  function  f-p(t)  whose 
Laplace  transform  is  <^)-p(s),  and  if  it  exists,  how  can  it 
be  determined  ?  The  answer  to  the  first  question  is  not 
always  positive.  In  general,  the  function  (p . j-(s)  must 
satisfy  some  restrictions  which  will  be  discussed  in  the 
next  subsection.  If  the  answer  to  ' the  first  question  is 
positive,  then  the  inversion  of  the  Laplace  transform  by  a 
complex  integral  gives  the  answer  to  the  second  question; 
That  is,  the  density  function  for  the  random  variable  T  is 
given  by  [Ref.  3] 

-I 

fT(t)  •  X[^T<S>] 

ot  +  i«o 

I  f  St  , 

- -TwrJ  e  •  9t(s)  ds- 

o<—  ioo 


19 


where  the  notation  JL  L9t(s)J  refers  to  the  inversion  of  the 
Laplace  transform  and  the  path  of  the  integral  above  is  the 
straight  line  Re  s=o<(o(>0)  parallel  to  the  imaginary  axis. 


D.  LAPLACE  INVERSION  AND  THE  THEOREM  OF  RESIDUES 

As  was  stated  in  last  subsection,  if  we  have  a  Laplace 
transform  <j6-p(s)  for  a  non-negative  random  variable  T,  then 
the  Laplace  inversion  T[<t>  y(s)j  gives  the  density  function 
f  T  (t). 

When  dealing  with  a  convolution  of  random  variables, 

that  is ,  T  =  T, +  . . .  +  Tn  ,  where  T( , . . . ,Tn  are  independent, 

the  convolution  density  function  is  given  by  [Ref.  3] 

-I 


f-r(t)  =  X(0t<s)] 


<rt  +  i  <*> 


■wr  J  e  •  0t<s)  ds> 


<*-  I  oo 


where  <^y(s)  is  the  Laplace  transform  of  the  sum 
(convolution)  T,  +  ...  +  T^  . 

To  obtain  the  convolution  density  function,  it  is  clear 
from  the  above  integral  that  we  have  to  evaluate  a  complex 
integral,  but  the  Theorem  of  Residues  reduces  the  evaluation 
of  a  complex  integral  along  closed  contours  to  passages  to 
the  limit  and  differentiations. 


Suppose  a  function  G(s)  is  holomorphic  in  a  finite 
domain  C,  the  term  holomorphic  meaning  that  the  function  is 


analytic  and  single-valued.  Then  the  integral 


1 

2  TT  i 


ds 


evaluated  along  it's  boundary,  vanishes.  If  the  function 
G(s)  is  not  holomorphic  in  C,  the  integral  may  not  be  zero. 
If  there  is  in  C  only  one  singularity  at  point  a,  the  value 
of  the  integral  is  called  the  residue  of  G(s)  at  the 


singular  point  a. 

The  well  known  Residues  Theorem  [Ref.  3]  is 


The  integral  of  an  analytic  function  along  a 
closed  contour  is  the  sum  of  the  residues  at 
the  singular  points  in  the  domain  enclosed  by 
the  contour,  multiplied  by  2  TT  i. 


Now,  let  us  define  the  analytic  function  G(s)  to  be 
G  (s )  =  e  S  ^  .  <j6T(s)  , 

where  <^>y(s)  is  the  Laplace  transform  of  the  sum  (convolu¬ 
tion)  T(  +  . . .  +  Tn  .  Then,  according  to  the  Theorem  of 
Residues,  the  convolution  density  will  be  obtained  as 
follows ; 


fT<‘)  - 

0<  +  '  00 

=  7777  /•  St  4>j(s)  ds 

<*-  i  eo 

at  +  i  oo 

■  tst  /G(s)  ds 

0<  —  I  oo 

m 

'  THfiJ  E/2"0  H  • 


where  m  is  number  of  singularities  (poles),  and  , 

i=l,...,m,  is  the  residue  of  the  analytic  function  G(s)  at 
the  singular  point  (pole)  a'L  . 

Therefore 

m 

'T(')  =  £  Ri  • 

i  - i 

that  is,  the  convolution  density  function  is  just  the  sum  of 
the  residue  of  the  analytic  function  G(s)  at  the  singular 
points . 

The  result  established  by  the  last  equation  would  not  be 
of  much  value  if  it  were  not  possible  to  compute  the  residue 
of  the  analytic  function  G(s)  at  the  singularities  directly 
and  without  evaluating  any  integral.  Thus,  let  us  now 
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This  implies  that  the  convolution  density  can  be  obtained  by 
limits  and  differentiations  which  are  easier  operations  to 
deal  with  than  the  integrations. 

Now,  we  are  ready  to  deal  with  the  three  different 
convolution  cases  which  have  been  mentioned  before.  we  will 
try  first  to  obtain  the  convolution  density  for  every  case. 
Then,  the  survival  functions  will  be  easily  obtained  by 
integration. 

E.  SURVIVAL  FUNCTION  OF  A  CONVOLUTION 

In  this  subsection,  we  will  use  the  results  obtained  in 
Subsection  D  to  find  the  survival  function  for  all  the 
possible  convolution  cases  which  have  been  mentioned  before. 
This  will  result  in  a  general  equation  for  the  survival 
function  of  any  convolution  of  independent  and  exponentially 
distributed  random  variables. 


T  =  T,  +  •••  +  Tn  , 

where  Tj  is  an  exponentially  distributed  random  variable 
with  a  failure  rate  X  i=l, . . . ,n  and  T, ,  ...  ,  Tn  are 
independent.  The  Laplace  transform  of  the  convolution  is 
given  by 


If  we  define  the  analytic  function 


G  ( s ) 


then,  according  to  the  Theorem  of  Residues,  the  convolution 
density  is  given  by 
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fT(t)  = 


Y  Ri  > 


where  n  is  the  number  of  singularities  (number  of  distinct 
failure  rates),  and  R  is  the  residue  of  the  analytic  func¬ 
tion  G(s)  at  the  singular  point  (pole)  -Xj  ,  i  =  1,...,  n. 

It  is  clear  from  the  Laplace  transform  function 
T(s)  that  there  are  n  singularities,  each  of  which  has 
an  order  (multiplicity)  one;  that  is,  nj  =  1  for  i  =  1,..., 
n.  Therefore,  the  residue  of  G(s)  at  the  singular 
point (pole)  -  X \  is  given  by 


R;  =  lim 


sdstoi 


=  lim 

s-^-X 


-J<-  *Xi 


)  G(s) 


By  substitution,  the  convolution  density  will  become 


f  T(t)  =  I- 


)  lim  , (s  +  X; )  e 

fr,  *— xf  ■ 


fj  (_V  ) 
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The  survival  function  of  the  system  will  be  obtained 
as  follows  ; 

Fy(t)  =  P  (  T  >  t  ) 
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2.  A  System  Having  n  Identical  Failure  Rates 


Consider  a  system  defined  by  the  shorthand  notation 
EXP (  \  )  +  ...  +  EXP (  \  )  (  n  times  ) 

The  life  time  of  the  system  is 

T  =  Tj  +  .  .  .  +  Tn  . 

where  T;  ,  i  =  l,...,n,  is  an  exponentially  distributed 
random  variable  with  a  constant  failure  rate  X  ,  and  T,  , . . . , 
Tn  are  independent.  The  Laplace  transform  of  the  convolu¬ 
tion  is 


0T(S) 


n  a o 

TT  f.-s* 


fT.  (t  )  dtf 


■vvv  w"."  v  "v  ‘trj  r 


T  -v  -  TV  V-v*y-i 


In  this  case  the  analytic  function  G(s)  has  one  sigular 
point  (pole)  of  order  (multiplicity)  n  at  the  point  s  =  -X  . 

According  to  the  Theorem  of  Residues,  the  residue  of 
the  analytic  function  G(s)  is  given  by 


*  fim 

s— 


J  ((S+\)«S))] 


Cn  -  0 1  dsn 

X" 


I 


=  Z//77  I  . 

S-,  Js 


--((S+X) 


(*  +  X)VJ 


I  </ 


rt-/ 


=  /  /  /77  .  .  . 

s— -U  (n-0!  Js 


S  f 


3=r(«  X  ) 


X 


n 


(n-  0 ! 


«*-'  -X£ 

t  e  ,  t>o. 


Therefore,  the  convolution  density  is  given  by 


f  T(  t )  =  R 


X 


Id 


(*  -  0  ! 


-Xt 

*  •  e 


-kt 


foT) 


t  •  e 


,  t  ^  0 , 
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which  is  the  density  function  of  the  Erlang  (n  ,  \  ) 
distribution. 


Now,  the  survival  function  is  given  by 
FT(t)  =  P  (  T  >  t  ) 


i )  du 


=  [f  («: 
t J 

eo  n  .  V 

r  x  -Xu 

=  J  — C± —  u  e 

.J  (r\-  Ol 


du 


Integrating  by  parts  yields 


Fy(t) 


A  (Xtf  -Xt 

L  (2-0'  e  '  *  >  0 

Z*l  '  1  ' 


which  is  the  well  known  formula  for  the  survival  function  of 
the  Erlang  (n  ,  X  )  distribution. 

3 .  A  Sytem  Having  Some  Similar  and  Some  Dissimilar 
Failure  Rate 

The  shorthand  notation  for  the  system  is 


EXP ( X,  )  +  ...  +  EXP (  X,  ) 

(n(  times) 

+  EXP(X2  )+...+  EXP  (  X2  ) 

(n 2  times ) 

+ 


+  EXP(\m)  +  ...  +  EXP(Xm)  (nm  times ) 

where  m  is  the  number  of  distinct  failure  rates,  and  n  is 
the  multiplicity  of  the  i-th  failure  rate  ,  i  =  l,...,m. 
The  Laplace  transform  of  the  convolution  is 


Then 


In  this  case,  the  system  has  m  distinct  failure  rates,  so 
that  the  analytic  function  G(s)  has  m  singularities. 
According  to  the  Theorem  of  Residues  the  convolution  density 
function  is  given  by 

m 

f  <t)  =  Y.  Ri  • 

is  I 
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where  R  is  the  residue  of  G(s)  at  the  singular  point  (pole) 

-  X  j  ,  i=  1 , . . .  ,  m. 

Since  n  j  is  the  multiplicity  for  the  failure  rate 
Xj,i  =  l,...,m,  it  follows  that  the  the  reside  of  the 
analytic  function  G(s)  at  the  singular  point  (pole)  -Xj  is 
given  by 


Rj  -  lim 
S 


lim 

S  -*  -X 


nj-i 

'  1  t 

.  ("i  -O' 

•  dsni-'V 

1 

n,-' 

^  (( 

(ni-0! 

ds^'-'V 

i 

J  ( 
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\'SrTT/  A  \J 
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,\l  (rtj  -  0!  <j. 


'  m  \ 

*  t  TT,  X  j 
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7 - TT  im  \  TFF'  W  +AJ) 
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ni*-/  *V-I 

Let  D  =  d/ds.  Then  D  =  d  /  ds  and 


.  (  TT  \?)  — lim  ,  e  *  *  RsA-f  1 


Let  us  define  the  function 


Then 


From  the  Binomial  Theorem  the  term  (D  +  t)  can  be  repre¬ 

sented  as 


n;  -/ 

(D  +  t) 

Therefore 
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Let  us  define  the  term 


C(i,k)  = 


r  d 


d  s 

K 


Hi (s1  ] 


•  -  Xi 


771 


-  TV, 


41-  TT(s  +  xJ)J 

J-K  T 


s  >; 


Then 


ffi  n;  ni'_l 


Ri  )  T. 


cu><) 


ki  c«i-K“Ol 


Since 


rn 


fT(t>  =  I 


1=1 


it  follows  that  the  convolution  density  is 


m  n.  m  .  n.-K-i 

t  \  \ v  co»  *)  + 


v*> 


j=l  Ks° 


<!  (nj-K-O! 


t  e 


This  equation  is  a  general  equation  for  the  density  function 
of  any  convolution  of  exponentially  distributed  random 
variables . 


Now,  the  survival  function  can  be  obtained 


follows ; 

FT(t)  =  p(  T  >  t  ) 


0 o 

]■ 


f  (u)  du 


("Mil  7^/“  «Xi“ 

I  =  I  Km0  •  v  1  '  t 


JlL 


m  r ..to  n;-l  it,-~  »  > 

i»i  K*o  ’At  Z»i  ^  ' 


Xit 


m  n.  m  n;-<  n;-* 


2-1 


cc>>>0  (Xjt)  —  t 


K,  (Z.,)!e 


,  t  >  o, 


i» l  K=o  2*» 


as 


This  is  a  general  equation  for  the  survival  function  of  any 
convolution  of  exponentially  distributed  random  variables. 


In  the  next  subsection,  we  will  give  three  examples 
of  the  application  of  the  equation.  These  examples  represent 
all  possiple  convolution  forms. 

F.  EXAMPLES  ON  THE  GENERAL  EQUATION 

1 .  A  Convolution  of  Dissimilar  Failure  Rates 

Using  the  general  equation  which  was  derived  in  the 
last  subsection,  this  example  illustrates  how  to  compute  the 
survival  function  for  a  system  defined  by  the  shorthand 
notation 

EXP(\,  )  +  EXP(X2)  +  EXP(X3) 

where 

a  i  ^  X  z  *  ^  i  ^  X3  *  X  2.  £  X  3  • 

Since  all  the  failure  rates  are  dissimilar,  it  follows  that 
the  total  number  of  distinct  failure  rates  is 

m  =  3  . 

The  multiplicity  of  the  first  distinct  failure  rate  X(  is 

n  j  =  1 . 

The  multiplicity  of  the  second  distinct  failure  rate  X2  is 
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The  multiplicity  of  the  third  distinct  failure  rate  \  is 

v 


nJ  =  1 


The  survival  function  of  the  system  is  given  by 


f  T  <t ,  . (fr\”)ZT  T* 


ial  K=OZ*l 


xr"  (*-»! 


=  suml  +  sum2  +  sum3 , 


where  suml  is  the  sum  at  i=l,  sum2  is  the  sum  at  i=2  and 
sum3  is  the  sum  at  i=3. 

At  i=l,  n,  =  1,  k  has  only  the  value  0,  and  z  has 
only  the  value  1.  The  value  of  the  derivative  term  C(1,0)  is 
given  by 


C(1,0) 


I 

•  ( s  t  X2 )  ( s  +  X3 ) 


s  = 


I 


(  X, )  (  X3-  X, ) 
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Therefore 


suml  = 


l-i  i-< 


X]  X3  Y,  X! 

K=o  2.-! 


C(i,*0  (X,t)  ^X|X 

\j~*  (2-0.' 


X  i  X  o  X- 


CU,o)  (X|  tY  £-X,t 
0 !  X,  0! 


X,  X,X3- - - - 

(  VX,)(X3-X,)X, 

Xz  X3  e-X,t 

(XZ-X,)(X3-X,) 


Similarly,  at  i=2,  =  1,  k  has  only  the  value  0, 

and  z  has  only  the  value  1.  The  derivative  term  C(2,0)  is 


and 


sum2 


1 

<2,0)  =  (\,-X*)(X,-Xz) 

X|  x3  e-X; 
( X|-  X2X  Xj-  X2) 


Also,  at  i=3,  nj  =  1,  k  has  only  the  value  0,  and  z 
has  only  the  value  1.  The  derivative  term  C(3,0)  is  given  by 


and 


k\ 

sum3  =  - - - e 

(  X,  -  X3)(  Xl  -  X3) 


.  x,t 


By  adding  suml,  sum2  and  sum3 ,  the  survival  function 


of  the  system  is 


F  (t)  =  suml  +  sum2  +  sum3 


X?  X’ 


X»  X2 


-X,t 


X,  X; 


_\3t 


(X2 - XD(X3-Xj)  (X^X^Xa-X^  Cx,- X3V)( X2-X^) 


,  t  >0 


2.  A  Convolution  of  Identical  Failure  Rates 


This  example  illustrates  how  to  compute  the  survival 
function  for  a  system  defined  by  shorthand  notation 


'(  X  )  +  EXP(  X  )  +  EXP(  X  ) 


This  convolution  has  only  one  distinct  failure  rate 
with  multiplicity  3.  Therefore  m=l  and  n=3.  The  survival 
function  of  the  system  can  be  obtained  as  follows; 


;=,  k=oZ=i  K X  ( z_,j • 
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U\ 
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r.-  ■. 
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r.  - 
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L 
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2  3-< 

K=0  Z*/ 


Z-t 


CO,  at)  (\t) 

k!  ■£-“  (z-o; 


Since  m  =  1,  and  n(  =  3,  it  follows  that  i  =  1  and  k 
0,1,2.  When  k  =  0  the  derivative  term  C(1,0)  is  given  by 


C(1,0) 


-*;■ 


-!  (s  +  \i) 


ds°  ^ 


—  X 


n(S.x/)nj  1 


Js=-X 


(S  +  0  )' 


,.-X 


Since 


Jsc 


#»♦  ] 
jf 


=  i, 


it  follows  that  the  k-th  derivative 


-  n : 


dK  /  dsKn(s+\j) 

ft' 


=  0  ,  k  1, 
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V  .  S.  , 


That  is,  C(l,k)  =  0,  for  k  =  1,2  .  Now,  the  survival  func¬ 
tion  for  the  system  will  become 


T<c>  =  X  I 


co,°)  (Xt)  -  X  t 

’  Xs  (*-')■' 


A  (Xt)"'  -Xt 

£(777)1  e 

2=1  X 


,  t  £  0, 


which  is  the  survival  function  equation  of  the  Erlang(  3,X) 


distribution . 


3.  A  Convolution  of  Similar  and  Dissimilar  Failure 
Rates 


This  example  illustrates  how  to  compute  the  survival 
function  for  a  system  having  some  similar  and  some  dissim¬ 
ilar  failure  rates,  using  the  general  equation. 


Suppose  we  have  a  system  defined  by  the  shorthand 


notation 


EXP (  X,  )  +  EXP(X,  )  +  EXPtXj.  ) 

In  this  case,  we  have  two  distinct  failure  rates  X|  and 
X2.  The  total  number  of  the  distinct  failure  rates  is 


m  =  2 
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The  first  failure  rate  has  the  multiplicity 


n 


i 

2. 


The  second  failure  rate 

n  2  = 


1  . 


has  the  multiplicity 


The  survival  function  equation  of  the  convolution  can  be 
given  as  follows; 


2  2 

ft(c)  -Or  xl)H  I  cc"*’ 


2-1 


/=!  K=0  Z*l 


*!  \?'K 


( X;  t) 

(z-0! 


-kit 


t  >  0. 


=  suml  +  sum2  , 


where  suml  is  the  sum  at  i=l,  and  sum^  is  the  sum  at  i=2. 

At  i=l,  k  takes  the  values  0  and  1.  We  will  evaluate 
the  two  derivatives  C(1,0)  and  C(l,l)  as  follows; 

2 


C(1,0) 


X 


I 


K=o  z«/  *  A. | 


Kt-'J  • 


At  i=2,  =1  and  k  takes  only  the  value  0  . 

we  have  to  evaluate  only  the  derivative  term  C(2,0) 
follows : 


III.  SURVIVAL  FUNCTION  AND  POLYNOMIAL  COEFFICIENTS 


In  section  2,  we  derived  a  general  equation  for  the 
survival  function  of  any  convolution  of  independent  and 
exponentially  distributed  random  variables.  This  section 
gives  an  alternative  formula  for  the  survival  function, 
using  the  idea  of  computing  the  coefficients  of  the  polyno¬ 
mials  that  accompany  the  exponential  terms  in  the  survival 

function  equation.  This  idea  was  explored  by  Sadan  Gursel 

[Ref.  4],  adding  one  distribution  at  a  time.  In  this  section 
we  will  derive  the  same  polynomial  coefficients  from  the 
general  equation  developed  in  the  previous  section. 

A.  DETERMINATION  OF  THE  POLYNOMIAL  COEFFICIENTS 

Consider  a  system  defined  by  the  shorthand  notation 
EXP(\|  )  +  ...  +  EXP(  X,  )  (n,  times) 

+  EXP  (X  ^  )  ♦  ...  -  EXP(X2)  (n2  times ) 

+  . 

+  EXP(Xjyj  )+  •••  +  EXP  ( X  jyj  )  (n777  times). 


This  notation  represents  the  convolution  of  n;  exponen¬ 


ts 


tial  random  variables,  where  there  are  n  identical  exponen¬ 
tial  random  variables  having  the  failure  rate 


X;  .  i  =  1,  •  •  •  » 


m  . 


We  expect  the  survival  function  of  the  convolution  to 
have  the  form 


-X,t 


■  Xjt 


-Xmt 


F  (t)  =  A,  (t )  .  e  +  A,(t).e  +  ...  +  Am(t).e  ,  t  =>  0 , 


where 


Aj(t)  -  a  ♦  a  t  ♦  a ,  ?  t 


n,-l 


I  z 


+  ai n.- \  c 


A2(t)  =  a2o  +  az\  t  +  a22  t  + 


r\z 


+  a2nr.  * 


2  *nT 1 

Am  (t)=  amo  +  ami  c  +  anu  c  +  •  •  *  +  am-n^-'t 

The  notation  Aj(t),  i=l,...,m,  represents  a  polynomial  of 
(nj-l)st  degree.  Therefore,  to  obtain  the  survival  function 


of  the  convolution,  we  have  to  find  all  the  coefficients 


For  example, 
notation 


the  convolution  represented  by  the  shorthand 


EX?  (  X  (  )  *  EXP  (  X ,  )  +  EXP(\,)  +  EXP(\2)  , 

where 

X  |  ^  X2  * 

has  a  survival  function  represented  by 


F  (t)  =  A ,  ( t )  .  e  X/  +  A2(t)  .  eX^^"  ,  t  >  0, 

where 

2. 

Al(t>  =  aIO  *  aM  c  *  al2  c  • 

and 

A 2 ( t )  =  ^  2.0 

It  is  important  to  note  that  the  number  of  exponential 
terms  is  equal  to  the  number  of  dissimilar  failure  rates, 
and  each  exponential  term  has  a  polynomial  coefficients  of 
degree  one  less  than  the  number  of  identical  random  vari¬ 
ables  having  the  corresponding  failure  rate. 

The  problem  of  computing  the  coefficients  will  be  solved 
if  we  derive  a  general  formula  for  the  constant  coefficients 
of  the  polynomial  coefficients  in  the  survival  function 
equation. 


Let  us  begin  with  the  general  equation  which  is 


Let  us  define 


G:(k,L)  = 
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C(i'K)  t 


K]  L  I 


Then, 
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+G j (n; -1,0) 
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*/- 1  (A--0-K 
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+G ; (nj-1,0) 


•G  (n: -2,1) 


+  G-  (0,2)+. . .+G; ( n  ■  -  3 , 2  ) 


+Gj (0 , n^ -1) 


Therefore 
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-X-t  n;  k  +  L  L 
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where 


"■•  ^cco*)  ,K+L  tL  • 
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Now,  it  easily  follows  that 


rn  i\\  L.  N  k+  L 

co,  <)  x. 


#11  Mi  '  • 

iL-<TJ,x.)E 


K  =0 


l  I  1 


represents  the  coefficient  of  t  for  the  i-th  polynomial  in 


the  survival  function  equation,  L  =  0,  ...  ,  n:-l, 


where 


C(i,k) 


This  formula  is  the  general  formula  for  the  coefficients  and 
it  can  be  considered  as  an  alternative  formula  for  the 
survival  function  of  any  convolution  of  exponential  random 
variables . 

It  is  clear  from  the  above  formula  and  from  the  general 
equation  in  Section  2  that  the  most  difficult  part  in  the 
reliability  computations  is  the  derivative  term  C(i,k), 
particularly  when  dealing  with  complex  convolutions  that 
have  a  long  stream  of  failure  rates  with  high  corresponding 
multiplicities.  We  will  give  the  the  practical  method  for 
computing  the  term  C(i,k)  in  Subsection  C. 


B .  EXAMPLE 

Consider  a  system  defined  by  the  shorthand  notation 
EXP  (  X ,  )  +  EXP(X2)  +  EXP  (  X  3  )  +  EXP(X3) 

Then,  the  survival  function  for  the  system  is  given  by 


F  (t) 


=  A, (t) 


e  1  + 


A3  (t ) 


A*(t).e  + 


The  first  distinct  failure  rate  \|  has  a  multiplicity  1, 
the  second  distinct  failure  rate  /\^  has  a  multiplicity  1, 
and  the  third  distict  failure  rate  A3  has  a  multiplicity  2. 
Therefore 


m  =  3  , 
n  |  =  1  , 
n^  =  1  , 

n3  =  2  . 

The  polynomial  Aj(t)  will  be  obtained  from 

"1-  1 

A,(t)  =  a|(J  ♦  a;|  t  ♦  ...  ♦  a ,  ,  t 

Substituting  n j  =  1  in  the  above  equation  yields 
Al (t)  =  a|0  • 

Since 
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it  follows  that 
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Now,  the  polynomial  A|(t)  is 
A,(t)  =  a  ,0 

2 

_ Xz.  X5 

"  (X2-X,)(X3-X,)z 

Similarly,  the  polynomial  A^(t)  is  given  by 
t  )  =  a  ^  0 
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X,  X 


l  ^3 


n3-l 


(X,  X Z )  (  X3  ”  X^) 

The  polynomial  A3(t)  will  be  obtained  from 

n 

A3(fc)  =  a30  +  a3l  t  +  •  •  •  +  *3^-1  fc 
Substituting  n^  =  2  yields 

A3^)  =  a30  +  a3\  t  ■ 

In  this  case  we  have  to  evaluate  the  two  coefficients 

a^0  and  a^j  .  Substituting  nj  =2,  m  =  3,  L  =  0  in  the 

general  equation  of  the  coefficients  yields 


"j  •  C(3,  K) 
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K.  Oi 
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X, 
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where 
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and 


3  _  -*j 

c(3,1)  '  [  4;  TJ(S  +  Xj' )  • 


^  (  X,  X?)  ( Xr"  Xa)  (  X,-  Xy)(X2~ \j) 


Thus 


X,  X?  X,  X2 X3 


X  |  X  2  X  3 


(X|  Xs)(  X2 "* X3)  (XrXa^Xjr)^)  (X(  Xj)(X2“X3) 


Substituting  L  =  1,  n^  =  2,  and  in  =  3  in  the  general  equation 
of  coefficients  yields 


( t  $ )  r  tt#  ^ 
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>,m  1  x. 


X|  X2.  ,0)  . 


Since 


C  ( 3 , 0  )  = 


•\b  irc^j''].. 
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it  follows  that 


'■'iw  « ■> 


(  X,-  \3)(  Xc  Xj) 


a3t  =  Xi  X2  Xr 


(Xr  Xj)(  X2-  X,) 


Finally,  the  survival  function  of  the  convolution  is 
given  by 


“X|t  , .  x  “X?t  .  ,.,~X3"t 


FT(t)  =  A  j  ( t )  .  e  1  +  A2(t).e  2  +  Aj(t)< 


a,0  e  ^  +  a7n  e  AZU  +  ( 


-X?t 


-X.t 


lZ  0 


a30  +a3l  O®  3 


X2Xj  eX,t 


X,  X3  e 


2-  -\zt 


( Xi~  X|)(  X3-  X,)  ( X|  X2)(X3  Xi) 


Xi  X2. 


X|  X2  X3 


(x~ X3)( Xz  X3)  ( X  - Xj) ( Xz  Xj) 


X|  X2  Xj 


XiX2Xjt  1  -X3t 


( X,-  X3)(  X2  X3 )  ( X,-  X3)(  X2  Xj) 


m 


56 


i  -  *  •  v  "  ^  m  m  j  ^  ‘  -  (  ,  ,  %  ^  ^  « 

*.*."**•«  m  *  ^  m  m  *  *  *  ,  *  ►  ■»  ^  k  "  ►  *  «  *  ,  *  n"'  r  ■  »  "  «  *  X  ’  •  "  »  *  *  *  »  ’  •>  */  *,*  *  .*  ,*  ,* 

1  /  ,^W.  A-’X--*  A..  A.”  a/a."  aT-VCi’  ^  •  -O-A."  /-  A  w  _r.  _ 


V  **A  V  A  *  A  * 


C.  THE  PRACTICAL  METHOD  FOR  COMPUTING  THE  DERIVATIVE  TERM 


C(I,K) 


One  of  the  major  difficulties  in  computing  the  survival 
function  either  by  the  general  equation,  or  by  the  general 
formula  of  the  coefficients,  is  the  derivative  term  C(i,k), 


where 


C(i,k)  = 


» 


This  term  represents  the  k-th  derivative  for  a  product  of 
m- 1  terms.  The  difficulty  in  comuting  C(i,k)  increases  with 
k.  Since  k  runs  from  zero  up  to  nj  -1,  i=l,...,m,  it  follows 
that  the  maximum  value  of  k  will  not  exceed  the  maximum 

I 

multiplicity  of  the  distinct  failure  rates  minus  one.  For 
example,  if  we  have  the  two  distinct  failure  rates  Xj  and 
Xz  where  X|  has  a  multiplicity  nj  and  ^as  a  multiplicity 
n^  ,  n(  >  n^  ,  then,  the  maximum  value  of  k  is  n, -  1. 

In  practice,  due  to  the  high  accuracy  (precision) 
required  in  the  computations,  we  will  limit  the  derivative 
term  C(i,k)  to  values  of  k  <9.  This  implies  that  the 
multiplicity  of  any  of  the  distinct  failure  rates  in  any 
convolution  must  not  exeed  10. 


The  following  is  one  of  the  solutions  which  will  be  used 
in  the  computer  program  for  computing  the  term  C(i,k)  (see 


Appendix.  A,  Appendix.  B). 

Suppose  we  have  a  convolution  which  has  m  distinct 
failure  rates  with  m  corresponding  multiplicities  .  We  know 
that 


Suppose  for  example  ,  we  want  to  compute  C(i,k)  at  i  =  1 , 
then  k  =  0,  ...  ,  nj-1,  and 


C(l,k) 


* 


il  ( _ 

dsK  Hst  XzY'1 


^  S  +  \m) 


s 


Our  problem  is  how  to  change  a  derivative  of  a  product  of 
terms  to  a  derivative  of  a  sum  of  terms,  since  the  later  is 
easier  in  computations. 

Let  us  define 


f 


TT  (s  + 


At  i=l  we  obtain 


-  n 


V  (s+\) 

i#' 


By  taking  the  natural  logarithim  on  both  sides,  we  obtain 


In  f[=  - 


X 

ln(s  *  Ay  )  + 


ln(s  + 


Now,  by  taking  the  first  derivative  on  both  sides  , 
obtain 


V  f.  = 


r  *2. 

1 

1  s*Xi+  + 

5  +  \m  J 

or,  equivalently 


/ 
f , 


I 

J*' 


"z 

1 

+  •  •  •  t 

L  s  +  Kz 

S  *  \m  J 

2  the  term 

"j 

s+  X; 

f, 
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f  «*v 


Then, 


m 

y  nJ 

s  +  X} 

j*1 

r  ryz 

|  1  -  _i_ 

nm  -j 

L  s  t  \z 

stL  -1 

Then 


f|=  -  A,  •  £, 


or,  equivalently 

df j /ds  =  -  A |  .  f j . 

When  i=l  and  k=l,  the  term  C(i,k)  is 

rn  -  fij 

ca,i>  •  [  -j-  TJ(s  +  \j)  ]  ^ 

l  As 

■  '  h  f.  ]  ...x, 


f,  =  -  A,  .  f,  , 


*  *  ■»  *  w  •*  ..  *  *  *» 

^  .  -  >.*  ,  -  .  ^ 
\-\v  •*  *  -* 


so  that 


C(l,l) 


We  can  proceed  in  a  similar  manner  to  compute  the 
derivatives  C ( 1 , 2 ) , C ( 1 , 3 ) , . . . , C ( 1 ,n  -1)  as  follows; 


^  / 
f ,  =  df,  /  ds 


d/ds 


-  A 


but 


so  that 


'll 


A,  f,  ♦  A,  f, 


2r 


=  (-D 


other 


where 


APPENDIX  A 


This  appendix  contains  a  computer  program  written  in 
Fortran  for  the  survival  function  of  any  convolution  of 
independent  and  exponentially  distributed  random  variables, 
using  the  general  equation 


rfi  n-  rtj- 

•*'  i«l  Ai  2=1 


i»l 


(Xjf) 

(z-i)! 


Xit  ts 

e  ,  t >o. 


A.  FUNCTIONS  AND  SUBROUTINES  USED  IN  THE  PROGRAM 

1.  FUNCTION  I FAC 

This  function  computes  any  factorial  required  in  the 

program. 

2.  SUBROUTINE  GAMMA 

This  subprogram  computes  the  incomplete  gamma  func¬ 
tion  with  the  parameters(  ,nj -k)  and  returns  the  result  in 
the  variable  G, 
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where 


3.  SUBROUTINE  DERIV 

This  subprogram  computes  the  k-th  derivative  C(i,k) 
for  k<9  and  returns  the  result  in  the  variable  C. 

B.  PROGRAM'S  LIMITATIONS 

(1)  The  multiplicity  of  any  of  the  distinct  failure  rates 
must  not  exceed  10 . 

(2)  The  computer  precision  limit. 

(3)  There  is  no  limitation  on  the  number  of  distinct  failure 


o 


rates . 


nononnnnnnnno 


section  Computer  program 


& 


$  JOB 


LAMDAI ....  I- th  distinct  failure  rate. 

LAMDA . Array  of  M  distinct  failure  rates. 

PRODCT ....  Product  of  all  the  failure  rates. 

FBAR . Survival  function  of  the  convolution. 

T . Time . 

C . Value  of  the  derivative  term  C(i,k). 

G . Incomplete  gamma  function. 

MULTI . Multiplicity  of  the  i-th  distinct 

failure  rate. 

MULTIP ... .Array  of  multiplicities  of  the  distinct 
failure  rates . 

M . #  of  distinct  failure  rates. 

I,J,K . Loop  index. 

DOUBLE  PRECISION  LAMDA(20),  LAMDAI , PRODCT , FBAR , 

*T,C,G 

INTEGER  I , J ,N ,M ,K ,MULTI , MULTIP (20 ) 

PRINT,  'Please  enter  the  time  of  duration' 

READ,  T 

PRINT,  'Please  enter  THE  #  OF  dissimilar  failure  rates’ 
READ,  M 

PRINT,  'Enter  the  values  of  dissimilar  failure  rates' 
READ,  ( LAMDA ( J )  ,  J=1,M) 

PRINT, 'Enter  corresponding  multiplicities 
READ,  (MULTIP (J)  ,J=1,M) 

PRINT , ' OUTPUT  ’ 

PRINT, ’======' 

WRITE ( 6 , 11)T 

FORMAT ('  TIME  =  '.F10.5) 

PRINT,'  FAILURE  RATE  MULTIPLICITY’ 

FBAR  =  0.0D0 
PRODCT  =  1.0D0 
DO  10  I  =  1,M 

WRITE ( 6 , 12 ) LAMDA ( I ) ,MULTIP(I) 

FORMAT (8X , F10 . 5 , 10X , 13 ) 

MULTI  =  MULTIP (I) 

PRODCT=PRODCT*LAMDA(I )**MULTI 
DO  20  J=l, MULTI 
K  =  J  -  1 
LAMDAI  =  LAMDA ( I ) 

CALL  GAMMA (MULTI, K,T, LAMDAI, G) 

CALL  DERIV ( I, M, LAMDA, MULTIP, K,C) 

FBAR=  FBAR+  C*G / ( I FAC ( K ) *LAMDAI** ( MULTI - K ) ) 

CONTINUE 

CONTINUE 

FBAR  =  FBAR  *  PRODCT 
WRITE ( 6 , 13 )T, FBAR 

FORMAT ( IX, 'P(T  >',F6.2,')  =  ’,F10.5) 

STOP 


MULTIPLICITY' 


••  •- .■ 

.  vl..' 
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END 

FUNCTION  IFAC(N) 

JFAC= 1 

IF(N.NE.O)  THEN  DO 
DO  30  J= 1 ,  N 
JFAC= JFAC*J 
30  CONTINUE 
END  IF 
IFAC= JFAC 
RETURN 
END 

SUBROUTINE  GAMMA (MULTI , K , T , LAMDAI , G ) 

DOUBLE  PRECISION  T,G, LAMDAI 
INTEGER  L,K,Z, MULTI 
G= .0D0 
L=MULTI-K 
DO  40  Z= 1 , L 

G  =  G + ( LAMDAI “T ) ** ( Z -  1 ) *DEXP ( - LAMDAI *T ) / 1 FAC ( Z -  1 ) 
40  CONTINUE 
RETURN 
END 

SUBROUTINE  DERIV (I ,M, LAMDA , MULTIP ,K , C ) 

DOUBLE  PRECISION  A , F , C , B ( 9 ) , LAMDA ( 10 ) 

INTEGER  I ,  J  ,  M ,  L ,  IK  , MULTIP ( 10 ) 

A= .0D0 
F= 1 . DO 
DO  50  J=I,M 
IF(I.NE.J)  THEN  DO 

F  =  F*  ( LAMDA ( J ) - LAMDA ( I ) ) ** ( - MULTIP ( J ) ) 

A=A+ (MULTIP (J ) / ( LAMDA (J ) -LAMDA(I ) ) ) 

END  IF 

50  CONTINUE 

IF (K . GT . I )  THEN  DO 
IK=K- I 

DO  60  L= 1 , IK 

B(L)= . 0D0 

DO  70  J=1,M 

IF(I.NE.J)  THEN  DO 

B (L) =B (L) + ( ( ( - 1 )**L)*MULTIP ( J )* 

* ( ( LAMDA ( J ) - LAMDA ( I ) ) ** ( - L -  1 ) ) ) *IFAC ( L ) 

END  IF 

70  CONTINUE 
60  CONTINUE 
END  IF 

IF(K.EQ.O)THEN  DO 
C  =  F 

ELSE  DO 

IF(K.EQ. 1)THEN  DO 
C=  (  -  1 )  *A*F 
ELSE  DO 

IF (K . EQ. 2 )THEN  DO 


^ .  ■  j  ■  i  ■»  >  ■■  f  -  » l!  -  '.  *  V  p  l  *J»  ■  l  ■  U  r<r  •< 


C= ( ( A**2 ) - B ( 1 ) ) *F 
ELSE  DO 

IF(K.EQ. 3)THEN  DO 

C= ( ( -  1 ) * ( A**3 )  +  3*A*B ( 1 ) - B ( 2 ) ) *F 

ELSE  DO 

IF(K.EQ.4)THEN  DO 

C= ( A**4 -  & *A** 2 *B ( 1 ) +4*A*B (2 )  +  3* 

*B ( 1 ) **2 - B ( 3 ) ) *F 
ELSE  DO 

IF (K . EQ . 5 )THEN  DO 

C= ( ( -  1 ) * ( A**5 )  + 10* (A**3 ) *B ( 1 ) -  10* ( A**2 ) * 

*b ( 2 ) - 15*A*B ( 1)**2+5*A*B ( 3 )  + 10*B ( 1 )*B (2 ) -B (4 ) )*F 


ELSE  DO 

IF (K . EQ . 6 )THEN  DO 
C  = ( A**  6-15  * ( A**4 ) *B ( 1 )  +  2  0* ( A**  3 ) * 

*B ( 2 ) -  15* (A**2 )*B ( 3 )  +  6*A*B (4 ) -B(5)  +  45 * 

* ( A**2 ) * ( B ( 1 ) **2 ) -  60 . DO*A*B ( 1 ) * 

*B ( 2 ) + 15*B ( 1 )*B ( 3 ) -  15 *B ( 1 )**3  +  10*B ( 2 ) **2 ) *F 
ELSE  DO 

IF(K. EQ. 7 )THEN  DO 

C  = ( ( -  1 ) * ( A** 7 )  +  2 i*a**5*B ( 1 ) -  3  5 *A**4*B ( 2 )  +  3  5 *A** 3  * 
*b ( 3 ) -  2 1*A**2*B (4 )  +  7*A*B ( 5 ) -B ( 6 ) - 105*A**3* 

*b ( 1 ) **2  +  2 10*A**2*B ( 1 ) *B ( 2 ) -  105 *A*B ( 1 ) *B ( 3 ) -  105* 

*B ( 1 ) **2  *B ( 2 ) -  7  0*a* 

*B ( 2 ) **2 +105  *A*B ( 1 ) **3  +  2 1*B ( 1 ) *B (4 )  +  3  5*B ( 2 ) * 
*B(3))*F 


ELSE  DO 

IF (K . EQ . 8 )THEN  DO 

C= ( A**8 - 28*A**6*B ( 1 ) + 56*A**5*B ( 2 ) - 70*A**4*B ( 3 ) + 56* 
A**3*B ( 4 ) - 28*A**2*B ( 5 ) +8*A*B ( 6 ) - B ( 7 ) + 2 10* 

A**4*B ( 1 ) **2 -  5  60*a**3*b ( 1 )* 

B ( 2 ) +420*A**2*B ( 1)*B (3 ) - 168*A*B ( 1 )*B (4 ) +28* 

B ( 1 ) *B ( 5 ) - 420*A**2*B ( 1 ) **3 + 840*A*B ( 1 ) **2*B ( 2 ) + 280* 
A**2*B ( 2 ) **2 -  2 10 *B ( 1 ) **2*B ( 3 ) - 280*A*B ( 2 )*B ( 3 ) - 280*B ( 1 ) 
B ( 2 ) **2  +  56  *B ( 2 ) *B (4 )  +  35*B ( 3 ) **2 + 105*B ( 1 ) **4 ) *F 
ELSE  DO 

IF (K . EQ . 9 )THEN  DO 

C= ( - A**9  +  36*A**7*B ( 1 ) - 84*A**6*B (2 )  +  126*A**5*B ( 3 ) -  126* 
*  A  * * 4  *  B ( 4 )  +  84*A**3*B ( 5 ) -  3  6*A**2*B ( 6 )  +  9*A*B ( 7 ) - B ( 8 ) - 
*3  7  8*A**5*B ( 1 ) **2 + 1230*A**4*B ( 1 ) *B ( 2 ) -  12  6  0*A**3*B ( 1 ) * 

*b ( 3 )  +  7  56*A**2*B ( 1 ) *B (4  ) -  2 52*A*B ( 1 ) *B ( 5 )  +  3  6*B ( 1 ) *B ( 6 ) ) 

C  =  C+ 1260 *A** 3  *B (1) **3-3780 *A**  2  *B ( 1 ) ** 2  * 

B ( 2 ) - 840*A**3*b ( 2 ) **2+ 18  90*a* 

B ( 1 ) **2*B ( 3 )  +  1260*A**2*B ( 2 ) *B ( 3 ) -  3  78*B ( 4 ) * 

B ( 1 ) **2 -404*A*B ( 2 ) *B ( 4 ) +  84*B ( 2 ) *B ( 5 ) - 945*a* 
b ( 1 ) **4 +  2  5  20*A*B ( 1 ) *B ( 2 ) **2  + 12  6  0*B ( 1 ) **3*B ( 2 ) -  3 15 *A* 
b ( 3 ) **2 - 1260*B ( 1 )*B (2 )*B ( 3 ) + 126*B ( 3 )*B (4 ) - 280*B ( 2 ) **3 
C=C*F 


END  IF 
END  IF 
END  IF 
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section  Examples  subsection  All  failure  rates  are  dissimil 
Shorthand  notation  : 

EXP (  .  3  )  +  EXP ( .32)  +  EXP ( . 4 )  +  EXP ( . 6  )  + EXP (.62) 

Output 

TIME  =  2.00000 

FAILURE  RATE  MULTIPLICITY 
0.30000  1 

0.32000  1 

0.40000  1 

0.60000  1 

0.62000  1 

P (T  >  2.00)  =  0.99817 
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subsection  All  failure  rates  are  identical 
Shorthand  notation  : 

EXP (2 . 1)  +  ...  +  EXP (2 . 1)  (10  times) 

Output 

TIME  =  3.20 

FAILURE  RATE  MULTIPLICITY 

2.10000  10 

P (T  >  3.20)  =  0.85772 


P  1  ■  ■  P  w  V  IF 


subsection  Some  similar  and  some  dissimilar  failure  rates 


Shorthand  notation  : 


♦  EXP (6.5)+ 

...  ♦  EXP (6. 5) 

(5  time 

EXP  (4 

•  5  )  +  ... 

+  EXP (4 . 5  ) 

(5 

t imes ) 

EXP  (3 

.5)*  ... 

+  EXP (3. 5) 

(5 

t imes ) 

EXP  (2 

•  5 ) +  ... 

♦  EXP (2. 5) 

(5 

t imes ) 

EXP  ( 1 

•  5  )  +  ... 

♦  EXP (1.5) 

(5 

times ) 

EXP  (  . 

50 ) ♦  . . . 

♦  EXP ( . 50 ) 

(5 

t imes ) 

Output 


TIME  =  8.00000 

FAILURE  RATE 

6.50000 

4.50000 

3.50000 

2.50000 

1.50000 
0.50000 

P (T  >  8.00) 


MULTIPLICITY 
5 
5 
5 
5 
5 
5 

0.99894 
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APPENDIX  B 


This  appendix  contains  an  alternative  computer  program 
written  in  Fortran  for  computing  the  survival  function  of 
any  convolution  of  independent  and  exponentially  distributed 
random  variables,  using  the  general  equation  for  the  coeffi¬ 
cients  of  the  polynomials  that  accompany  the  exponential 
terms  in  the  survival  function  equation 


F(t)  = 


,  t  >0  , 


where 


a 


it 


m 


"j 


jfi 


I 


K-0 


CU,k) 
K!  Li 


K+  L 

X{ 


A.  FUNCTIONS  AND  SUBROUTINES  USED  IN  THE  PROGRAM 

1 .  Function  IFAC 

This  function  computes  any  factorial  required  in  the 


program . 


2.  Subroutine  COEFF 


This  subprogram  computes  the  coeffient  of  t  for  the 
i-th  polynomial  in  the  survival  function  equation  and 
returns  the  result  in  the  variable  CO, 


where 


m  "i-'-L 

»■  (TjM)Z 

jfi 


CC « >  *) 
<i  L ! 


K  +  L 


X, 


3 .  SUBROUTINE  DERIV 


This  subprogram  computes  the  k-th  derivative  C(i,k) 
for  k<9  and  returns  the  result  in  the  variable  C. 


B.  PROGRAM'S  LIMITATIONS  : 

(1)  The  multiplicity  for  any  of  the  distinct  failure  rates 
must  not  exceed  10. 

(2)  The  computations  must  be  within  the  computer  precision 
limit . 


(3)  No  limitation  on  the  number  of  distinct  failure  rates. 


onnnnononoooon 


section  Computer  program 


$  JOB 

LAMDAI . . . . i- th  distinct  failure  rate. 


LAMDA . Array  of  M  distinct  failure  rates. 

MULTI . Multiplicity  of  the  i-th  distinct 

failure  rate. 


MULTIP ... .Array  of  multiplicities  of  the  distinct 
failure  rates . 

PRODCT ....  Product  of  the  distinct  failure  rates. 


FBAR . Survival  function  of  the  convolution. 

T . Time . 

C . Value  of  the  derivative  term  C(i,k). 

CO . Coeffetiont  of  the  polynomials. 

M . #  of  distinct  failure  rates. 

I,J,K . Loop  index. 


DOUBLE  PRECISION  LAMDA (20) ,  LAMDAI , PRODCT , FBAR , 

*T , C , CO 

INTEGER  I , J , N ,M , K , MULTI , MULTIP (20 ) 

PRINT,  'Please  enter  the  time  of  duration' 

READ,  T 

PRINT,  'Please  Enter  #  of  dissimilar  failure  rates' 
READ,  M 

PRINT,  'Enter  the  values  of  dissimilar  failure  rates' 
READ,  ( LAMDA ( J )  ,  J=1,M) 

PRINT, 'Enter  corresponding  multiplicities  ' 

READ,  (MULTIP (J)  ,J=1,M) 

PRINT , ' OUTPUT  ' 

PRINT, ’======' 

WRITE (6, 11 )T 

11  FORMAT ('  TIME  =  ',F10.5) 

PRINT,'  FAILURE  RATE  MULTIPLICITY' 

FBAR  =  0.0D0 
DO  10  I  =  1,M 

10  WRITE (6,12) LAMDA ( I ) , MULTIP ( I ) 

12  FORMAT (8X , F10 . 5 , 10X , 13 ) 

DO  20  I  =  1 ,M 

MULTI  =  MULTIP (I ) 

WRITE (6 ,13)1 

13  FORMAT (IX, 'POLYNOMIALS  COEFFICIENT  FOR  LAMDA’, 12,’  is 
DO  30  LL  =  1, MULTI 

L  =  LL  -  1 

CALL  COEFF ( I, L,M, MULTIP, LAMDA, CO) 

WRITE ( 6 , 14 ) I , L , CO 

14  F0RMAT(3X, 'a' ,11,11, '  =  ’.F25.12) 

FBAR= FBAR+ T**L*C0*DEXP ( - LAMDA ( I ) *T ) 

30  CONTINUE 
20  CONTINUE 

WRITE (6, 15 )T, FBAR 

15  FORMAT ( / , IX , ' P (T  >',F6.2,')  =  ’.F20.12) 
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s  % 


STOP 

END 

FUNCTION  IFAC(N) 

JFAC=  1 

IF(N.NE.O)  THEN  DO 

DO  40  J=  1 ,  N 

JFAC= JFAC*J 

CONTINUE 

END  IF 

IFAC= JFAC 

RETURN 

END 

SUBROUTINE  COEFF ( I ,  L ,  M , MULTIP , LAMDA , CO ) 

DOUBLE  PRECISION  C , CO , LAMDA ( 10 ) , PRODCT 
INTEGER  I , L , M , K , KK , LL , MULTIP ( 10 ) 

PRODCT  =  1.0D0 

DO  50  J  =  1,M 

IF(I.NE.J)  THEN  DO 

PRODCT = PRODCT *LAMDA  ( J )  —MULTIP  (  J  ) 

END  IF 

CONTINUE 

CO  =  O.ODO 

LL  =  MULTIP ( I )  -  L 

DO  60  KK  =  1 , LL 

K  =  KK  -  1 

CALL  DERIV ( I, M, LAMDA, MULTIP, K,C) 

CO  =  CO  +  C*PRODCT*LAMDA(I)**(K+L)/IFAC(K)/IFAC(L) 

CONTINUE 

RETURN 

END 

SUBROUTINE  DERIV ( I , M , LAMDA , MULTIP , K , C ) 

DOUBLE  PRECISION  A ,  F  ,  C  ,  B  ( 9  ) , LAMDA ( 10 ) 

INTEGER  I , J ,M,L , IK , MULTIP ( 10 ) 

A= .0D0 
F= 1 . DO 
DO  70  J= 1 ,M 
IF(I.NE.J)  THEN  DO 

F  =  F* ( LAMDA ( J ) - LAMDA (I) ) ** ( - MULTIP ( J ) ) 

A=A+ (MULTIP ( J) / ( LAMDA ( J ) - LAMDA (I ) ) ) 

END  IF 
CONTINUE 

IF(K.GT.l)  THEN  DO 
IK=K- 1 

DO  80  L= 1 , IK 
B (L) = . 0D0 
DO  90  J= 1 ,M 
IF(I.NE.J)  THEN  DO 
B(L)=B(L)+(((- 1)**L)*MULTIP( J)* 

*( (LAMDA(J)-LAMDA(I) )**(-L-l) ) )*IFAC(L) 

END  IF 
CONTINUE 


/  .*  .•  •  -•  .• 
•V.V- 


.N  .*• 


V-V-y  \ 


•-'.'••Vv'-V. 


•  r  *  V  .* 


CONTINUE 
END  IF 

IF(K.EQ.O)THEN  DO 
C=F 

ELSE  DO 

IF(K. EQ. 1)THEN  DO 
C= ( - 1)*A*F 
ELSE  DO 

IF(K.EQ.2)THEN  DO 
C= ( (A**2) -B ( 1) )*F 
ELSE  DO 

IF(K.EQ.3)THEN  DO 

C= ( ( - 1 ) * ( A**3 )  +  3*A*B ( 1 ) - B ( 2 ) ) *F 

ELSE  DO 

IF(K.EQ.4)THEN  DO 

C= (A**4- 6*A**2*B ( 1 ) +4*A*B ( 2 )  +  3* 

*B(1)**2-B(3))*F 
ELSE  DO 

IF(K. EQ. 5 )THEN  DO 

C= ( ( - 1 ) * ( A**5 )  + 10* ( A**3 )*B ( 1 ) -  10* ( A**2 ) * 

*B ( 2 ) -  15  *A*B ( 1 ) **2  +  5  *A*B ( 3 ) + 10*B ( 1 ) *B ( 2 ) -B (4 ) ) *F 
ELSE  DO 

IF(K.EQ.6)THEN  DO 

C= (A** 6  - 15* (A**4 ) *B ( 1 )  +  20* (A**3 ) * 

*B ( 2 ) -  15* ( A**2 ) *B ( 3 )  +  6*A*B ( 4 ) -B ( 5 )  +  45* 

* ( A** 2 ) * ( B ( 1 ) **2 ) -  6  0 . DO*A*B ( 1 ) * 

*B  (  2  )  + 15 *B  ( 1 )  *B  (  3  )  -  15 *B  ( 1 ) **3  +  10*B  ( 2  )  **2  )  *F 
ELSE  DO 

IF(K.EQ. 7 )THEN  DO 

C= ( ( -  1 ) * ( A**7 )  +  2 1*A**5*B ( 1 ) -  3  5*A**4*B ( 2 )  +  3  5*A**3* 

*B ( 3 ) - 21*A**2*B (4 )  +  7*A*B ( 5 ) -B ( 6 ) - 105*A**3* 

*B ( 1 ) **2  +  2 10*A**2*B ( 1 ) *B ( 2 ) - 105*A*B ( 1 ) *B ( 3 ) -  105* 

*B ( 1 ) **2*B ( 2 ) - 70*A* 

*B ( 2 ) **2 +105  *A*B ( 1 ) **  3  +  2 1 *B ( 1 ) *B ( 4 )  +  3  5  *B ( 2 ) * 

*B ( 3 ) )*F 
ELSE  DO 

IF (K. EQ. 8 )THEN  DO 

C= (A**8 - 28*A**6*B ( 1 )  +  5  6*A**5*B ( 2 ) -  7  0*A**4*B ( 3 )  +  56* 
*A**3*B ( 4 ) -  2  8  *A**2*B ( 5 )  +  8  *A*B (6)-B(7)  +  210* 

*A**4*B ( 1 )**2- 560*A**3*B ( 1 )* 

*B ( 2 ) +420*A**2*B ( 1 ) *B ( 3 ) - 168*A*B ( 1 ) *B ( 4 ) +28* 

*B ( 1 ) *B ( 5 ) - 420*A**2*B ( 1 ) **3  +  840*A*B ( 1 ) **2*B ( 2 )  +  280* 
*A**2*B ( 2 ) **2 -  2 10*B ( 1 ) **2*B ( 3 ) - 280*A*B ( 2 ) *B ( 3 ) - 280*B ( 1 ) 
*B ( 2 ) **2  +  5  6  *B ( 2 ) *B ( 4 )  +  3  5*B ( 3 ) **2  + 10  5  *B ( 1 ) **4 ) *F 
ELSE  DO 

IF(K . EQ. 9 )THEN  DO 

C= ( - A**9  +  36*A**7*B ( 1 ) - 84*A**6*B ( 2 )  +  126*A**5*B ( 3 ) -  126* 
*A**4*B ( 4 )  +  84*A**3*B ( 5 ) -  3  6  *A**2*B ( 6 )  +  9*A*B ( 7 ) - B ( 8 ) - 
*  3  7  8  * A**  5  *B ( 1 ) ** 2  + 1 2  3  0  *A**4 *B ( 1 ) *B ( 2 ) -  1 2  6  0 *A**  3  *B ( 1 ) * 
*B ( 3 )  +  756*A**2*B ( 1 )*B (4 ) - 252*A*B ( 1 )*B (5 )  +  36*B ( 1 )*B (6 ) ) 
C=C+1260 *A** 3 *B ( 1 ) ** 3  -  3 7 8 0 *A**2*B ( 1 ) ** 2* 


*B ( 2 ) - 840*A**3*B ( 2 ) **2  +  18  90*A* 

*B ( 1 ) **2*B ( 3 )  + 12  6  0*A**2*B ( 2 ) *B ( 3 ) -  3  78*B ( 4 ) * 

*B ( 1 ) **2 - 404*A*B ( 2 ) *B ( 4 )  +  84*B ( 2 ) *B ( 5 ) - 945*A* 

*B ( 1 ) **4  +  25  2  0*A*B ( 1 ) *B ( 2 ) **2  + 12  6  0 *B ( 1 ) ** 3  *B ( 2 ) -  3 1 5 * A* 
*B(3)**2-1260*B(1)*B(2)*B(3)+126*B(3)*B(4)-280*B(2)**3 
C=C"F 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
RETURN 
END 

$ ENTRY 


section  Examples  subsection  All  failure  rates  are  dissimilar 
Shorthand  notation  : 

EXP ( . 3 )  +  EXP ( . 32 ) +EXP ( . 4 ) +EXP ( . 6 ) +EXP ( . 62 ) 

Output 


TIME  =  2.00000 

FAILURE  RATE  MULTIPLICITY 
0.30000  1 

0.32000  1 

0.40000  1 

0.60000  1 

0.62000  1 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  1  IS  : 

alO  =  248.000000000000 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  2  IS  : 

a20  =  -332.142857142857 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  3  IS  : 

a30  =  101.454545454545 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  4  IS  : 

a40  =  -70.857142857143 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  5  IS  : 
a50  =  54.545454545454 

P (T  >  2.00)  =  0.998171698511 
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subsection  All  failure  rates  are  identical  : 
Shorthand  notation  : 

EXP (2 . 1 )  +  ...  +  EXP (2.1)  (10  times 

Output 


TIME  =  3.20000 

FAILURE  RATE  MULTIPLICITY 

2.10000  10 

POLYNOMIAL  COEFFICIENT  FOR  LAMDA  1  ARE  : 


alO 

= 

1.000000000000 

all 

= 

2.100000000000 

al2 

= 

2.205000000000 

al3 

= 

1.543500000000 

al4 

= 

0.810337500000 

al5 

= 

0.340341750000 

al6 

= 

0.119119612500 

al7 

= 

0.035735883750 

al8 

- 

0.009380669484 

al9 

= 

0.002188822880 

(T  > 

3.20)  = 

0.857717853017 

subsection  Some  similar  and  some  dissimilar  failure  rates 


Shorthand  notation  : 


EXP (6.5)+  . 

,  .  + 

EXP (6 .5) 

(5 

times ) 

EXP (4.5)+  . 

,  .  + 

EXP (4 . 5 ) 

(5 

times ) 

EXP (3.5)+  . 

,  .  + 

EXP (3. 5) 

(5 

times ) 

EXP (2.5)+  . 

,  .  + 

EXP (2. 5) 

(5 

times ) 

EXP (1.5)+  . 

,  .  + 

EXP (1.5) 

(5 

times ) 

EXP (.50)+  . 

,  .  + 

EXP ( . 50) 

(5 

times ) 

Output 


TIME  =  8.00000 


FAILURE  RATE 

6.50000 

4.50000 

3.50000 

2.50000 

1.50000 
0.50000 


MULTIPLICITY 

5 

5 

5 

5 

5 

5 


POLYNOMIALS  COEFFICIENTS  FOR  LAMDA  1  ARE  : 
alO  =  -0.033655633163 

all  =  -0.015984638484 

al2  =  -0.002968083834 


a!3 


-0.000255685577 


al4 


-0.000008633539 


POLYNOMIAL  COEFFICIENTS  FOR  LAMDA  2  ARE 
a20  =  78824.407500664440 

a21  =  28606.965700048750 

a22  =  4242.082627828356 

a23  =  308.298280647481 

a24  =  9.469913057431 

POLYNOMIAL  COEFFICIENTS  FOR  LAMDA  3  ARE 
a30  =  -1939257.342698666000 

a3 1  =  -833405.255639648600 

a32  =  -199021.555234616700 

a33  =  -18294.927183023870 

a34  =  -1641.852439502142 

POLYNOMIAL  COEFFICIENTS  FOR  LAMDA  4  ARE 
a40  =  2591869.114320258000 

a41  =  -558729.118591437400 

a42  =  340359.641443397200 

a43  =  -14083.275902489360 

a44  =  4142.139971320400 

POLYNOMIAL  COEFFICIENTS  FOR  LAMDA  5  ARE 
a50  =  -736153.955891001200 

a51  =  310180.360064797400 

a52  =  -58447.530039485250 


a53 


5362.158719218832 


a54  =  -297.897706623269 

POLYNOMIAL  COEFFICIENTS  FOR  LAMDA  6  ARE 


a60  = 

4718.810424379766 

a61  = 

-1542.762826749752 

a62  = 

203.749271368373 

a63  = 

-12.977288263887 

a64  = 

0.350737520646 

P(T  >  8.00) 


0.998940885595 
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